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Abstract 

Applications  utilizing  smart  materials  are  rapidly  increasing  and  include  high  speed  milling  and  hybrid  motor 
design.  Such  application  utilize  magnetostrictive  transducers  operating  in  hysteretic  and  nonlinear  regimes.  To 
achieve  the  high  performance  capabilities  of  these  transducers,  models  and  control  laws  must  accommodate  the 
nonlinear  dynamics  in  a  manner  which  is  robust  and  facilitates  real-time  implementation.  To  this  end,  the 
models  and  control  algorithms  must  utilize  known  physics  to  the  highest  degree  possible,  be  low  order,  and  be 
sufficiently  robust  to  operate  under  realistic  conditions.  In  this  paper  we  consider  the  robust  control  of  a  smart 
structure  with  disturbances  due  to  inherent  hysteresis  and  sensor  noise.  We  demonstrate  the  techniques  on  a 
magnetostrictive  transducers  but  they  are  sufficiently  general  to  be  utilized  on  several  commonly  used  smart 
materials.  The  performance  of  the  control  strategies  are  illustrated  through  numerical  examples. 
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1.  Introduction 

Smart  materials  are  increasingly  utilized  for  a  variety  of  applications,  including  tunable  lenses,  modular 
antennas,  atomic  force  microscopes,  and  high  speed  milling.  Materials  such  as  piezoceramics  (PZT)  and  magne- 
tostrictives  are  ideal  for  such  applications  because  of  their  set  point  accuracy  and  high  bandwidth  capabilities. 
However,  these  materials  exhibit  hysteresis  and  constitutive  nonlinearities  which  must  be  accommodated  to 
achieve  the  high  performance  demands  of  these  applications.  At  low  frequencies  and  moderate  drive  levels, 
these  effects  can  often  be  mitigated  through  feedback  loops.  At  high  drive  levels  or  high  frequencies,  however, 
the  hysteresis  and  nonlinear  dynamics  must  be  incorporated  into  models  and  subsequent  control  designs.  One 
method  of  attenuating  the  nonlinear  hysteretic  behavior  of  these  materials  is  to  develop  an  inverse  model.  This 
method  requires  a  hysteresis  model  which  allows  an  inverse  that  can  be  calculated  in  real-time.  This  type  of 
inverse  filtering  has  been  applied  previously  employing  a  domain  wall  hysteresis  model  [4,  5]  as  well  as  Preisach 
models  [10]  (for  a  description  of  the  domain  wall  hysteresis  model,  see  [1,  2]).  In  this  paper  we  develop  an  inverse 
compensation  method  which  utilizes  a  free  energy  hysteresis  model.  In  a  realistic  application,  the  effects  of 
hysteresis  are  not  the  only  disturbances  to  be  overcome.  There  is  a  certain  amount  of  noise  associated  with  any 
sensing  device.  This  noise  will  be  fed  into  the  controller  and  can  greatly  decrease  the  closed  loop  performance.  In 
this  paper,  we  consider  the  development  of  a  robust  feedback  control  law  which  is  able  to  reject  the  disturbances 
due  to  sensor  noise  and  the  hysteretic  behavior. 

To  illustrate,  we  will  consider  the  prototypical  magnetostrictive  actuator  depicted  in  Figure  1.  Input  stresses 
and  displacements  are  provided  by  the  Terfenol-D  rod  in  response  to  fields  generated  by  the  surrounding  solenoid. 
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Figure  1.  Terfenol-D  transducer. 


The  structure  is  encased  in  a  permanent  magnet  which  provides  a  biasing  field  to  obtain  bi-directional  strains 
and  a  mechanism  for  flux  shaping  in  the  transducer.  As  detailed  in  [2],  such  actuators  have  the  capability 
of  generating  broadband,  high  force,  responses.  However,  they  also  exhibit  varying  degrees  of  hysteresis  and 
nonlinearities  in  the  relation  between  the  input  field  H  and  the  magnetization  M  in  the  Terfenol-D  rod.  We 
concentrate  on  a  magnetostrictive  transducer  in  this  paper  but  note  that  the  hysteresis  model  and  control  design 
presented  are  sufficiently  general  to  permit  direct  extension  to  analogous  piezoelectric  and  ferroelectric  models 
as  developed  in  [7,  8,  9]. 


2.  Free  Energy  Hysteresis  Model 


In  this  section  we  shall  summarize  the  free  energy  hysteresis  model  presented  in  [6] .  This  model  is  based  on  the 
quantification  of  energy  required  to  reorient  dipoles  in  combination  with  stochastic  homogenization  techniques 
to  accommodate  variations  in  coercive  and  effective  fields.  For  this  development,  we  assume  fixed  temperatures. 
The  model  also  ignores  losses  due  to  eddy  currents  and  therefore  should  be  employed  for  low  frequency  drive 
levels.  A  further  assumption  is  that  the  spring  washer  in  the  transducer  provides  enough  prestress  to  dominate 
the  crystal  anisotropies  (see  Figure  1). 

We  first  consider  the  Helmholtz  energy  which  incorporates  the  internal  energy  due  to  the  interaction  of  dipole 
moments.  Under  the  assumption  that  dipoles  orient  either  with  the  applied  field  or  diametrically  opposite,  it  is 
illustrated  in  [6]  that  a  reasonable  form  of  the  Helmholtz  energy  is 
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As  depicted  in  Figure  2,  Mr  and  Mj  respectively  denote  the  point  at  which  the  minimum  of  "0  occurs  and  the 
inflection  point. 

In  the  presence  of  an  applied  field  H,  the  Gibbs  energy  is  given  by 


G  =  'ip-  HM. 


(2) 


Note  that  the  magnetostatic  energy  is  e  =  hqHM  where  fig  denotes  the  magnetic  permeability.  The  Gibbs  free 
energy  given  in  (2)  can  be  interpreted  as  incorporating  /io  into  ip  for  simplicity. 

For  a  homogeneous  material  with  the  effective  field  taken  as  the  applied  field,  the  local  average  magnetization 
can  be  expressed  as 


M  =  x+{M+)  +  x_{M_) 


(3) 


¥(M)=  G(0,M) 


Figure  2.  (a)  Helmholtz  energy  ip  and  Gibbs  energy  for  increasing  field  H.  (b)  Local  magnetization  M  for 
increasing  field  H. 


where  x+  and  X-  denote  the  fraction  of  moments  having  the  positive  or  negative  orientation,  respectively,  and 
(M_|_)  and  {M-)  are  the  expected  values  of  the  magnetization  for  moments  in  the  positive  or  negative  orientation 
respectively.  As  detailed  in  [6],  (M_|_)  is  quantified  by 

dM 

e-G(H.M)v/kT 

Here  g-GiH,M)V/kT  i-epj-esentg  probability  of  obtaining  an  energy  level  G  and  the  denominator  of  (4)  assures 
that  the  probability  of  integrating  over  all  possible  magnetization  values  is  1.  Here  V  is  the  lattice  volume,  k  is 
Boltzman’s  constant  and  T  is  the  temperature.  The  relation  for  (Af_)  is  given  similarly  as 


(M_) 


J  —OO 

g-G(H,M)V/kT 

J  — CXD 


(5) 


The  moment  fraction  are  quantified  by  the  evolution  equations 


x+  =  -P+-X+  +  P-+X-  (6) 

X- = —p — ^x-+p+-x+  (7) 

where  p^ _ and  p _ are  the  likelihoods  of  switching  form  positive  to  negative  orientation  and  switching  from 

negative  to  positive  orientation,  respectively.  They  are  given  as 
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where  m  is  the  mass  of  lattice  volume  V.  The  relation  between  the  applied  field  H  and  the  magnetization 
M  exhibits  both  hysteresis  and  nonlinear  transition  because  the  local  magnetization  (3)  is  probabilistic.  The 
steepness  of  the  transition  depend  on  the  ratio  of  GV  to  kT. 

We  will  now  describe  an  alternative  formulation  of  the  magnetization  which  we  will  employ  numerically.  We 


Figure  3.  Varying  field  levels  for  the  limiting  case  local  average  magnetization  M. 


can  analyze  the  equilibrium  behavior  after  the  moments  have  switched  by  examining  the  equilibrium  condition 
=  0  which  yields  H  =  Therefore  the  local  model  predicts  a  linear  relationship  between  H  and  M 

with  a  slope  of  K  We  can  employ  asymptotic  relations  to  simplify  the  magnetization  and  consider  the  jumps  in 
magnetization  to  occur  instantaneously.  As  detailed  in  [6],  the  local  magnetization  in  this  limiting  case  is  given 

as  _ 

(  [MiH;H,,Om  ,T(t)=0 

[M{H;  =  I  f  ~  ^  0  and  iJ(max  T{t))  =  -H^  (10) 

[  ^  +  Mr  ,  T{t)  yf  0  and  iJ(max  r(t))  =  He 
where  the  transition  points  are  specified  as 

T{t)  =  {<  e  (0,  Tf]  I  Hit)  =  -He  or  H{t)  =  He}  (11) 

and  the  initial  moment  orientation  is  given  as 

r  f  -  Mr  ,  i/(0)  <  -He 

\M{H-He,mO)=  l  ^  .-He<H{0)<He  (12) 

[  f  +  Mr  ,  H{0)  >  He  . 

Figure  3  depicts  the  progression  of  local  average  magnetization  points  for  varying  applied  field  levels. 

The  relation  (10)  is  derived  under  the  assumption  that  the  lattice  structure  is  homogeneous  and  hence  the 
domain  structure  is  homogeneous.  This  implies  that  the  free  energy  profile  for  the  different  regions  of  the 
Terfenol-D  are  identical.  This  assumption  is  overly  simple  since  it  ignores  material  defects,  polycrystallinity  and 
nonuniformities  in  the  crystals.  Also,  the  model  has  assumed  that  the  effective  field  at  the  domain  level  is  the 
applied  field.  To  extend  the  local  model,  stochastic  distributions  are  incorporated  to  develop  a  bulk  magnetization 
model  for  a  nonhomogeneous  Terfenol-D  sample  with  nonconstant  effective  fields  (see  [6]  for  details). 

We  shall  begin  by  summarizing  the  inclusion  of  a  distribution  of  free  energy  profile  to  accommodate  material 
nonhomogeneities  (see  Figure  4).  The  nonhomogeneities  are  included  by  assuming  the  parameters  Mr  and  Mj, 
or  equivalently  He  =  r]{MR  —  Mj),  are  normally  distributed  with  mean  He-  The  total  magnetization  can  then 
be  given  as 

poo 

M{H)=  M{H-,Hc,0f{He)  dHe  (13) 

Jo 

with  density 

f{He)  =  (14) 

The  parameters  Ci  and  b  are  positive  and  M  is  given  by  (3)  or  (10). 

The  second  extension  entails  the  incorporation  of  variations  in  the  effective  field.  The  effective  field  is  assumed 
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Figure  4.  (a)  Free  energies  for  nonhomogeneous  domain  structure,  (b)  Resulting  variations  in  the  local  coercive 
field  He- 

to  be  normally  distributed  with  mean  TL.  Therefore,  for  a  fixed  coercive  field,  the  magnetization  can  be  given  as 

/OO  _ 

M{H;Hc,0C2e-^^-‘^^"^^dndHc.  (15) 
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Combining  coercive  field  and  effective  field  distributions  for  the  full  magnetization  of  a  nonhomogeneous  poly¬ 
crystalline  sample  of  Terfenol-D  with  variable  effective  fields  yields  the  magnetization  relation 

POO  POO  _ 

M{H)  =  C  /  (16) 
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For  numerical  implementation  of  the  free  energy  hysteresis  model,  the  distributions  for  the  effective  field  and 
coercive  field  are  evaluated  using  a  composite  quadrature  rule.  This  allows  the  hysteresis  model  to  be  computed 
algebraically  and  hence  at  low  cost.  Because  of  the  exponential  decay  of  the  distributions,  they  can  be  truncated 
to  speed  computation.  The  model  is  sufficiently  simple  to  allow  the  possibility  of  real  time  implementation  and 
has  a  relatively  low  number  of  material  dependent  parameters. 


3.  Inverse  Hysteresis  Model 

To  construct  an  inverse  compensator,  the  monotonicity  of  the  hysteresis  model  is  exploited.  To  determine 
the  magnetic  field  required  to  create  a  desired  magnetization,  the  hysteresis  model  is  advanced  until  the  desired 
magnetization  is  surpassed.  Then  the  magnetic  field  is  computed  by  a  linear  interpolation  between  the  last  two 
points.  The  computational  speed  of  the  inverse  compensator  depends  on  the  size  of  the  step  taken  in  advancing 
the  hysteresis  model.  Larger  steps  will  increase  the  speed  while  decreasing  the  accuracy  of  the  inverse  compen¬ 
sator.  To  facilitate  real  time  control,  the  control  method  must  be  designed  to  reject  a  significant  amount  of  error 
in  the  linearization  of  the  hysteresis.  The  linearization  error  for  an  input  signal  with  a  frequency  of  1  Hz  and  a 
step  size  of  AiJ  =  1  employed  in  the  inverse  model  is  plotted  in  Figure  5a.  While  the  error  appears  quite  large 
the  relative  error,  illustrated  in  Figure  5b,  is  reasonable.  The  control  design  must  be  able  to  reject  this  error  to 
the  input  of  the  plant. 


4.  Full  Actuator  Model 


We  have  described  the  model  for  the  magnetization  in  the  Terfenol-D  sample  due  to  an  applied  field.  We 
must  next  quantify  the  strains,  forces,  and  displacements  generated  by  these  changes  in  magnetization  for  the 


prototypical  transducer  under  consideration  (see  Figure  1).  The  model  for  the  transducer  is  based  on  the  model 
given  in  [2]. 

One  end  of  the  rod  {x  =  0)  is  assumed  to  be  fixed  while  the  other  end  {x  =  L)  is  constrained  by  a  damped 
oscillator  and  has  an  attached  point  mass  (see  Figure  6).  The  internal  damping  coefficient,  density,  Young’s 
Modulus  and  point  mass  are  given  by  cd,  P,  E,  and  respectively.  The  damping  spring  has  stiffness  and 
Kelvin- Voigt  damping  coefficient  cl- 

Under  the  assumption  of  linear  elasticity,  Kevin- Voigt  damping  and  small  displacements,  the  stress  at  any 
point  X,  0  <  a;  <  L,  is  given  as 

a{t,x)  =  +  -  E—M{t,xf  (17) 


where  w{t,x)  is  the  displacement  and  Ag  and  Mg  respectively  denote  the  saturation  free  strain  and  saturation 
magnetization  of  the  Terfenol-D.  Balancing  forces  yields 
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where  the  resultant  is  specified  by 
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To  obtain  appropriate  boundary  conditions,  we  first  note  zc(t,  0)  =  0.  Balancing  forces  at  x  =  L  yields 

L)  =  -kLw{t,  L)  -  L)  -  L). 
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Figure  6.  Rod  model. 


Figure  7.  Robust  control  diagram. 


Initial  conditions  are  taken  to  be  ry(0,a;)  =  0  and  (0,a;)  =  0. 

ox 

The  model  can  be  implemented  utilizing  the  weak  form  and  a  Galerkin  finite  element  approximation.  To 
further  simplify  our  system,  we  assert  that  the  magnetization  of  the  Terfenol-D  rod  can  be  taken  as  uniform 
over  the  length  of  the  rod.  This  is  reasonable  since,  in  present  actuator  design,  flux  shaping  via  the  surrounding 
magnet  can  be  used  to  minimize  end  effects  in  the  rod.  Hence  motion  of  the  tip  of  the  Terfenol-D  rod  can  be 
modeled  by  an  ordinary  differential  equation.  This  observation  can  be  supported  by  the  fact  that,  since  each 
section  of  the  rod  reacts  identically  to  the  uniform  magnetic  field,  the  elements  act  uniformly  and  therefore  the 
transducer  dynamics  may  be  modeled  as  a  damped  spring-mass.  The  ODE  model  we  employ  for  the  displacement 
of  the  tip  of  the  Terfenol-D  rod  due  to  the  magnetization  is  given  as 

mx  +  kx  +  cx  =  ujM‘^{t). 

x(0)  =  xq  x(0)  =  xo-  (21) 

Here  the  scalars  m,k,c,and  oj  are  determined  by  fitting  the  model  (21)  to  the  Galerkin  approximation  of  (18)  or 
to  data  from  the  physical  device. 

Next  we  shall  linearize  the  magnetostrictive  relationship  about  a  biasing  magnetization  level  of  Ms/2.  This 
bias  can  be  achieved  my  means  of  the  permanent  magnet  (see  Figure  1).  The  dynamics  of  the  Terfenol-D 
transducer  can  be  now  represented  by 

X  +  kx  +  cx  =  ujM  {t) 

x(0)  =  Xo  x(0)  =  Xo-  (22) 

where  k  =  k/m  =  7.8899  x  10^,  c  =  c/m  =  6.4251  x  10^  and  w  =  1.3724  x  10“^.  Note  that  the  hysteresis  inherent 
in  the  Terfenol-D  rod  is  still  present  in  (22)  by  means  of  the  hysteretic  relationship  between  H  and  M. 


5.  Robust  Control  Design 

In  this  section,  we  present  a  control  design  which  incorporates  the  presence  of  external  disturbances  and 
minimizes  their  effects  with  respect  to  a  given  norm.  We  shall  consider  the  H2  norm  defined  by 

1  k°° 

ni^  =  — y  trace[T*(ja;)r(ja;)]d^.  (23) 

The  H2  norm  is  chosen  because  it  provides  a  measure  of  energy  in  the  system.  The  inverse  compensator  devel¬ 
oped  in  Section  3  will  be  employed  to  attenuate  the  hysteresis  inherent  in  the  Terfenol-D  transducer.  Figure  7 
shows  the  block  diagram  of  the  system  to  be  controlled.  In  the  diagram  P  represents  the  transducer  model  given 
by  the  differential  equation  (22).  The  signal  to  be  tracked  is  denoted  by  r  while  the  position  of  the  tip  of  the 
Terfenol-D  rod  is  denoted  by  y.  The  signal  d  represents  the  error  in  the  linearization  of  the  input  by  the  inverse 
filter  (see  Figure  8).  The  signals  n  and  s  represent  noise  in  the  sensor  the  measurements  of  y.  For  the  purposes 
of  this  paper,  we  assume  sensor  noise  with  a  frequency  of  60  Hz.  This  noise  is  represented  by  s.  The  signal  n  rep¬ 
resents  higher  frequency  noise  which  could  be  attributed  to  the  sensing  device,  other  external  disturbances  or,  in  a 


Figure  8.  Inverse  compensation. 


numerical  simulation,  numerical  noise.  We  have  chosen  to  separate  these  two  signals  as  we  wish  to  weight  them 
independently.  The  output  signals  e  and  u  denote  the  weighted  tracking  error  and  weight  output  of  the  controller 
K,  respectively.  Here  Wd,  We,  Wr,  Ws  and  Wn  are  weighting  functions  which  are  chosen  to  maximize  the 
performance  of  the  controller  utilizing  a  priori  knowledge  of  the  characteristics  of  the  signals. 

We  can  now  represent  the  transfer  functions  of  the  open  loop  system.  The  maps  from  the  inputs,  r,  d,  s  and 
n,  to  the  outputs  u,  e  and  v  are  given  as 

V  =  Wr[r]  -  iP[Wd[d]  +  u]  +  WM  +  WJs])  =  W,[r]  -  P[Wd[d]]  -  WM  -  Wjs]  -  PM 


e  =  We  [t>] 
u  =  Wu[u]  . 

Denoting  the  transfer  function  matrix  from  the  inputs  to  the  outputs  as  G,  we  then  have 


G  = 


WeWr  -WePWd  -WeWn  -WeWs  -WeP 

0  0  0  0 

Wr  -PWd  -Wn  -VFs  -P 


(24) 


(25) 


The  system  can  be  represented  as  a  linear  fractional  transformation  as  shown  in  Figure  9.  The  goal  is  to  design 
a  controller  K  so  the  P2  norm  of  transfer  function  from  the  inputs  to  the  outputs  of  our  closed  loop  system, 
denoted  T,  is  minimal. 


The  open  loop  system  can  be  partitioned  as 


A 

Pi 

P2 

G{s)  = 

Cl 

0 

Di2 

.  C2 

D21 

0 

(26) 


where 


Gii  = 


A 

Pi  ' 

Cl 

0 

Gi2  = 


A 

P2 

Cl 

P12 

G21  = 


A 

Pi 

C2 

D21 

G22  = 


A 

P2  ' 

C2 

0 

(27) 


represent  the  transfer  functions  from  w  to  z,  u  to  z,  w  to  v,  and  u  to  v,  respectively  (see  Figure  9).  We  make 
the  following  assumptions  on  our  system: 

1.  {A,  Bi)  is  controllable  and  (Gi,  A)  is  observable. 

2.  {A,  B2)  is  stabilizable  and  (G2,  A)  is  detectable. 


4. 


12P12  >  0  and 

A 

-juji 

P2 

Cl 

P12 

A 

-jujI 

Pi 

C2 

P21 

has  full  column  rank  for  all  oj. 


5. 


has  full  row  rank  for  all  uj. 
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Figure  9.  Linear  fractional  transformation  (LFT)  representation  of  transducer. 


5.1  Weighting  Functions 

The  majority  of  effort  in  designing  a  robust  controller  lies  in  the  choice  of  the  weighting  functions.  A  discussion 
on  the  choice  of  these  function  can  be  found  in  [3,  11].  The  frequency  of  the  noise  s  present  can  be  accurately 
determined  for  the  device  measuring  the  position  of  the  tip  of  the  Terfenol-D  rod.  The  weight  Ws  is  chosen 
to  heavily  weight  frequencies  in  a  bandwidth  around  our  expected  noise  frequency  and  attach  low  weights  to 
frequencies  above  and  below  our  specified  bandwidth.  For  the  simulation,  a  noise  frequency  of  60  Hz  is  added  to 
the  system  and  Ws  is  taken  to  be  a  ninth-order  pass-band  filter  with  a  bandwidth  of  10  Hz  centered  at  60  Hz.  A 
higher  order  filter  can  be  used  if  a  steeper  rolloff  outside  of  the  frequency  band  is  desired  although  we  have  chosen 
to  limit  the  number  of  states  in  the  state  space  representation.  The  frequency  response  of  Ws  is  illustrated  in 
Figure  10a. 

To  weight  high  frequency  noise,  Wn  is  taken  to  be  a  high  pass  filter.  The  sixth-order  filter,  plotted  in 
Figure  10b,  is  employed  in  the  numerical  simulations.  The  weight  Wr  is  taken  as  a  sixth-order  pass-band  filter. 
Since  the  frequency  of  the  reference  signal  r  is  known,  we  can  design  Wr  to  have  a  narrow  bandwidth  which  is 
taken  as  1  Hz  centered  at  1  Hz.  The  frequency  response  of  Wr  is  shown  in  Figure  11. 

To  determine  the  weighting  function  W^,  a  signal  with  the  same  frequency  as  the  reference  signal  was  filtered 
by  the  inverse  compensator  and  then  fed  into  the  free  energy  hysteresis  model.  The  output  is  defined  to  be  the 
desired  control  signal  plus  the  disturbance  d  (see  Figure  8).  The  Fourier  transform  of  the  disturbance  is  analyzed 
to  determine  Wd-  Figure  12a  illustrates  the  power  spectrum  of  the  disturbance.  From  Figure  12a,  it  is  clear  that 
Wd  should  weight  frequencies  between  100  and  350  Hz  heavily.  Therefore  Wd  is  taken  as  a  pass-band  filter  with 
a  bandwidth  of  250  centered  at  225  Hz.  The  frequency  response  of  Wd  is  shown  in  Figure  12b. 

The  weighting  function  on  the  error  signal  was  taken  to  be  We  =  — —  with  -je  =  3.2  x  10®  and  Ee  =  1  x  10“®. 


(a)  (b) 

Figure  10.  Frequency  response  of  (a)  Ws  and  (b)  Wn- 


An  integrator  was  chosen  to  prevent  the  error  from  achieving  steady  state  at  a  nonzero  value  and  the  pole  was 
shifted  slightly  off  zero  to  ensure  the  controller  designed  is  realizable.  We  take  the  weighting  function  on  the 
controller  output  to  be  Wu  =  5  x  10“®.  Since  we  do  not  experience  any  problems  with  saturation,  or  other  such 
effects,  we  minimally  weight  u  to  focus  the  controller  on  tracking  and  disturbance  rejection. 


5.2  H2  Optimal  Control  Design 


Employing  the  notation  defined  in  (27),  the  design  of  the  optimal  H2  controller  incorporates  two  Hamiltonian 
matrices 


and 


Ho  = 


Ji  = 


A-  BoR^^DloCi  -B2H[^BI 

-Cl{I  -  Di2Ri^Dlo)Ci  -(A  -  B2R^^DloCi) 


{A  -  B2R^^DloCiY 


2-^2  ^2 


(28) 


(29) 


-Bi{I  -  D*oiR2^D2i)Bt  -{A  -  B2R^^DtoCi) 

where  Ri  =  DI2D12  >  0  and  R2  =  D21D21  >  0. 

These  Hamiltonian  matrices  give  rise  to  two  Riccati  equations 

{A  -  B2R^^Dl2Ci)*X2  +  X2(A  -  B2Ri^D*2Ci)  +  X2(-B2Ri^B;)X2  -  Cl{I  -  Di2R^^Dl2)Ci  (30) 


and 


(a) 


(b) 


Figure  12.  (a)  Power  spectrum  of  d,  and  (b)  frequency  response  of  Wd- 


{A  -  B2R^^Dl^Ci)Y2  +  Y2{A  -  B2R^^Dl2Ciy  +  Y2{-C;R^^C2)Y2  -  Bi{I  -  21) Bl.  (31) 

Now  that  the  transducer  is  represented  and  the  Riccati  equations  have  been  defined,  the  following  theorem  from 
[11],  guaranteeing  an  H2  optimal  controller,  can  be  employed. 


Theorem  1:  There  exists  a  unique  controller  which  minimizes  the  H2  norm  of  the  closed  loop  system  if: 

1.  H2  Gdom(Ric)  and  X2  :=Ric(i/2)  >  0 

2.  J2  €dom(Ric)  and  Y2  :=Ric(J2)  >  0 


The  H2  optimal  controller  is  given  by 


where 


A2 

—L2 

F2 

0 

(32) 


A2  =  A  +  B2F2  +  L2C2,  F2  =  -Ri\BIX2  +  Dl^Ci),  L2  =  -{Y2CI+BiDy)R, 


2^- 


(33) 


5.3  Numerical  Results 

The  numerical  results  for  the  H2  robust  control  law  are  presented  in  this  section.  The  results  were  computed 
using  the  noise  signal  s  with  a  magnitude  of  1  x  10“®  and  a  frequency  of  60  Hz.  The  inverse  compensator  was 
computed  using  4  quadrature  points  and  20  divisions  for  the  coercive  field  and  40  divisions  for  the  effective  field. 
A  step  size  of  AiJ  =  1  was  employed  for  the  hysteresis  inverse  model.  These  parameters  ensure  a  simulation  of 
the  control  process  which  has  the  potential  for  real-time  implementation.  The  tracking  performance  with  the 
presence  of  sensor  noise  and  the  disturbance  is  illustrated  in  Figure  13.  Figure  13a  illustrates  the  controller’s 
ability  to  track  the  reference  signal  and  reject  the  noise  and  disturbance  signals.  A  tracking  error  less  that  2 
microns  is  achieved  after  a  short  period,  as  illustrated  in  Figure  13b,  and  the  noise  signals  and  effects  of  the 
inverse  compensation  for  the  hysteresis  are  effectively  attenuated. 


(a)  (b) 

Figure  13.  H2  (a)  tracking  performance  and  (b)  tracking  error. 


6.  Conclusions 


This  paper  has  outlined  an  H2  robust  control  design  for  a  smart  system.  Although  the  methods  were  illus¬ 
trated  on  a  magnetostrictive  transducer,  they  can  be  applied  to  systems  utilizing  piezoceramic  or  shape  memory 
alloys  as  well.  It  was  shown  that  an  H2  control  design  was  capable  of  maintaining  accurate  tracking  while 
rejecting  sensor  noise  and  a  disturbance  due  to  the  inexact  inverse  filter.  The  same  methodology  can  be  used 
to  include  other  disturbances  to  the  system  if  needed.  Details  were  provided  regarding  the  choice  of  weighting 
functions,  or  filters,  since  this  step  in  the  robust  control  formulation  is  crucial  to  achieving  the  balance  between 
tracking  requirements  and  disturbance  rejection  needed  for  a  particular  application. 
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